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We present measurements of the branching fraction B and longitudinal polarization fraction /l 
for ai(1260)+ ai(1260)" decays, with ai(1260)* ■k'-k+'k^. The data sample, collected 

with the BABAR detector at the SLAC National Accelerator Laboratory, represents 465 x 10® 
produced B'B pairs. We measure ^ ai(1260)+ ai(1260)- ) x [B(ai(1260)+ 7r-7r+7r+) f = 

(11.8 ± 2.6 ± 1.6) X 10"® and /l = 0.31 ± 0.22 ± 0.10, where the first uncertainty is statistical and 
the second systematic. The decay mode is measured with a significance of 5.0 standard deviations 
including systematic uncertainties. 



PACS numbers; 13.25.Hw, 12.15.Hh, ll.30.Er 
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Charmless B decays to final states involving two axial- 
vector mesons (AA) have received considerable theoret- 
ical attention in the last few years P, 0]- Using QCD 
factorization, the branching fractions of several B AA 
decay modes have been calculated. Predictions for the 
branching fraction of the B° ai(1260)+ ai(1260)" de- 
cay mode vary between 37.4 x 10~^ [l| and 6.4 x 10^^ 
0. Branching fractions at this level should be observ- 
able with the BaBAR data sample, which can be used 
to discriminate between the predictions. The predicted 
value of the longitudinal polarization fraction is 0.64 
The only available experimental information on this 
B decay mode is the branching fraction upper limit (UL) 
of 2.8 X 10"'^ at 90% confidence level (CL) measured by 
CLEO 0. 

The measured value /l ~ 0.5 in penguin-dominated 
B (j)K* decays [3| is in contrast with naive standard 
model (SM) calculations predicting a dominant longitu- 
dinal polarization (/^ ^ \) 'm B decays to vector-vector 



(VV) final states. The naive SM expectation is confirmed 

pp and B+ ^ ujp+ [6\ 



1 is found in vector-tensor 
while /l ~ 0.5 is found in 
(see Ref [sl for further dis- 



in the tree-dominated B 
decays. A value of Jl 
B (j)K*{U30) decays 
B LoK*{U30) decays 
cussion). 

The small value of /i observed in i? — > (pK* decays 
has stimulated theoretical effort, such as the introduc- 
tion of non-factorizable terms and penguin-annihilation 
amp litudes Other explanations invoke new physics 
[lOj . Measurement of /l in B'^ af Gi [ii[ decays will 
provide additional information. 

We present the first measurements of the branching 



fraction and polarization in B*^ 



decays, with 



|12| . We do not separate the P-wave 
(7r7r)p and the S-wave (7r7r)cr components in the ai Stt 
decay; a systematic uncertainty is estimated due to the 
difference in the selection efficiencies . Due to the lim- 
ited number of signal events expected in the data sample, 
we do not perform a full angular analysis. Using helicity 
formalism, and after integration over the azimuthal an- 
gle between the decay planes of the two oi mesons, the 
predicted angular distribution dT/d cos 9 is: 



1 dT 



Pdcos 61 



i/T(l + cos2 0), (1) 



where fx = II and 9 is the angle between the normal 
to the decay plane of the three pions of one ai and the 
fiight direction of the other ai, both calculated in the rest 
frame of the first ai. 

The results presented here are based on data collected 
with the BaBAR detector lj| at the PEP-II asymmetric- 



energy e+e collider [15| located at the SLAC National 



Accelerator Laboratory. The analysis uses an integrated 
luminosity of 423.0 fb~^, corresponding to (465±5) x 10^ 
BB pairs, recorded at the T(4S') resonance at a center- 
of-mass energy of y/s — 10.58 GeV. An additional 



43.9 fb-\ taken about 40 MeV below this energy (off- 
resonance data), is used for the study of qq contimmm 
background (e+e~ qq , with q = u,d, s, c). 

Charged particles are detected, and their momenta 
measured, by a combination of a vertex tracker (SVT) 
consisting of five layers of double-sided silicon microstrip 
detectors, and a 40-layer central drift chamber (DCH), 
both operating in the 1.5 T magnetic field of a super- 
conducting solenoid. The tracking system covers 92% of 
the solid angle in the center-of-mass frame. We identify 
photons and electrons using a CsI(Tl) electromagnetic 
calorimeter (EMC). Further charged-particle identifica- 
tion is provided by the specific energy loss (dE/dx) in 
the tracking devices and by an internally reflecting ring- 
imaging Cherenkov detector (DIRC) covering the central 
region. A K/ n separation of better than four standard 
deviations is achieved for momenta below 3 GeV/c, de- 
creasing to 2.5 a at the highest momenta in the B decay 
final states. A more detailed description of the recon- 
struction of charged tracks in BaBAR can be found else- 
where [lij . 

Monte Carlo (MC) simulations of the signal decay 
mode, continuum, BB backgrounds and detector re- 
sponse fr^ are used to establish the event selection crite- 
ria. The MC signal events are simulated as B^ decays to 
a'^ with ai — > p(770)7r. The oi meson parameters in 
the simulation are: mass uiq = 1230 MeV/c^ and width 
To = 400 MeV/c^ 

We reconstruct the decay of into three charged pi- 
ons. Two pion candidates are combined to form a p'^ 
candidate. Candidates with an invariant mass between 
0.51 and 1.10 GeV/c^are combined with a third pion to 
form an ai candidate. The ai candidate is required to 
have a mass between 0.87 and 1.75 GeV/c^. We impose 
several particle identification requirements to ensure the 
identity of the signal pions. We also require the prob- 
ability of the B vertex fit to be greater than 0.01 and the 
number of charged tracks in the event to be greater or 
equal to seven. 

A B meson candidate is kinematically charac- 
terized by the energy-substituted mass toes = 

(s/2 + pg • pb)'^/Eq — p| and energy difference AE = 
Eg — y^/2, where the subscripts and B refer to the ini- 
tial T(4S') and the B candidate in the laboratory frame, 
respectively, and the asterisk denotes the T(4S') frame. 
The resolutions in mES and AE are about 3.0 MeV/c^ 
and 20 MeV, respectively. We require candidates to sat- 
isfy 5.27 < mES < 5.29 GeV/c^ and -90 < A£; < 70 
MeV. 

Background arises primarily from random track com- 
binations in continuum events. We reduce this back- 
ground by using the angle 9t between the thrust axis 
of the B candidate and the thrust axis of the rest of the 
event (ROE), evaluated in the T(4S') rest frame. The 
distribution of | cos 9t \ is sharply peaked near 1 for com- 
binations drawn from jet-like continuum events and is 
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nearly uniform for BB events; for this reason, we require 
I cos 6*1! < 0.65. 

Background can also arise from BB events, especially 
events containing a charmed meson (these are mostly 
events with five pions and a mis-identified kaon in the 
final state). The charmed background includes peaking 
modes, with structures in TOes and AE that mimic sig- 
nal events, and non-peaking "generic" modes. To sup- 
press the charm background, we reconstruct D and D* 
mesons. Events are vetoed if they contain D or D* candi- 
dates with reconstructed masses within 20 MeV/c^ (win- 
dow size of about ±2cr) of the nominal charmed meson 
masses [l8| . 

The mean number of B candidates per event is 2.9. If 
an event has multiple B candidates, we select the can- 
didate with the highest B vertex probability. From 
MC simulation, we find that this algorithm selects the 
correct candidate 90% of the time in signal events while 
inducing negligible bias. 

Using MC simulation of signal events with longitudinal 
(transverse) polarization, signal events are divided in two 
categories: correctly reconstructed signal (CR), where all 
candidate particles come from the correct signal S", and 
self-cross feed (SCF) signal, where candidate particles 
are exchanged with a ROE particle. The fraction of SCF 
candidates is 31.8 ± 3.2 (19.4 ± 1.9)%. 

We determine the number of signal events (the signal 
yield) from an unbinned extended maximum-likelihood 
(ML) fit. The seven input observables are AE, tubs, a 
Fisher discriminant J-^ jl6| . the two Oi masses and the 
two H — \cos9\. The Fisher discriminant combines 
four variables calculated in the T(45') frame: the abso- 
lute values of the cosines of the angles with respect to 
the beam axis of the B momentum and the thrust axis 
of the B decay products, and the zeroth and second angu- 
lar Legendre moments Lo,2 of the momentum flow about 
the B thrust axis. The Legendre moments are defined 
by ifc = I^mPm |cos6'm|'', where 9m is the angle with 
respect to the B thrust axis of a track or neutral cluster 
m, Pm is its momentum, and the sum includes the ROE 
particles only. 

There are five hypotheses in the likelihood model: sig- 
nal, continuum, and three BB components, which take 
into account charmless, generic charm and peaking charm 
backgrounds. The likelihood function is: 



N 



5 

E 



(2) 



where N is the number of input events, rij is the number 
of events for hypothesis j and Vj{'x.i) is the corresponding 
probability density function (PDF), evaluated with the 
observables of the ith event. Since correlations among 
the observables are small (< 10%), we take each V as the 
product of the PDFs for the separate variables. 



The signal includes both CR and SCF signal compo- 
nents with the SCF fraction fixed in the fit to the value 
estimated from MC simulation. Both CR and SCF sig- 
nals are used to measure the branching fraction and po- 
larization. The PDF of the signal takes the form: 



p 

^ stg 



- fL{l-gl^'n'PcR,L + fL9f^''^VsCKL (3) 
+ /t (l — 9t"^) '^CR.T + frgT^^'PsCF.T 

where (?£ (St^^) is the fraction of SCF in longitudinal 
(transverse) polarized signal events and Vcr,l, 'Pscf,l 
[VcR.T, Vscf,t) are the signal PDFs of CR and SCF 
signal components for longitudinal (transverse) polariza- 
tion. 

We determine the PDF parameters from Monte Carlo 
simulation for the signal and BB backgrounds and from 
off-resonance data for the continuum background. 

We parameterize toes and AE using a Gaussian 
function with exponential tails for the CR sig- 

nal and charmless components, and using polynomi- 
als for all other components, except for the TOes 
distribution for continuum events which is described 
by the ARGUS empirical phase space function 21 1 
— exp [—^(1 — a;^)] , where x = 2m-E,s/ \fs and ^ 
is a parameter. The ai mass is described by a relativis- 
tic Breit-Wigner function for the CR signal component, 
an asymmetric Gaussian plus a linear polynomial for the 
SCF signal component, and polynomials for the remain- 
ing components. The Fisher variable is parametrized 
with an asymmetric Gaussian plus a linear polynomial 
in all cases. The Ti variables are parametrized with a 
Gaussian plus a linear polynomial for the charm peak- 
ing component and with a polynomial in all other cases. 
The parameters left free in the fit are the signal, contin- 
uum, and three BB component yields, and /l. We also 
float some of the parameters of the continuum PDFs: the 
three parameters of the asymmetric Gaussian part of T , 
and one parameter each for the 7i, the ai masses and 
AE. 

Large data samples of B decays to charmed final states 
[B^ — > D*~a^), which have similar topology to the sig- 
nal, are used to verify the simulated resolutions in toes 
and AE. Where the data samples reveal differences from 
the Monte Carlo we shift or scale the resolution function 
used in the likelihood fits. Any bias in the fit, which arises 
mainly from neglecting the small correlations among the 
discriminating observables, is determined from a large set 
of simulated experiments for which the continuum back- 
ground is generated from the PDFs, and into which we 
have embedded the expected number of BB background, 
signal and SCF events chosen randomly from fully simu- 
lated Monte Carlo samples. 

The fit results are presented in Table HI The detection 
efficiencies are calculated as the ratio of the number of 
signal MC events passing all the cuts to the total num- 
ber generated. We compute the branching fraction by 
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TABLE I: Fitted signal yield and yield bias (in events), bias 
on /l, detection efficiencies el and er for events with longi- 
tudinal and transversal polarization, respectively, significance 
S (including systematic uncertainties), measured branching 
fraction B and fraction of longitudinal polarization Jl with 
statistical and systematic uncertainties. 



Signal yield 

Signal yield bias 

fibias 

eL (%) 

er (%) 

S{a) 

B [xlQ-^) 
h 



545 ± 118 
+14 
-0.06 
9.0 
10.0 
5.0 

11.8 ± 2.6 ± 1.6 
0.31 ±0.22 ±0.10 



subtracting the fit bias from the measured yield, and di- 
viding the result by the number of produced BB pairs 
times the product of the daughter branching fractions 
and the detection efficiency. We assume that the branch- 
ing fractions of the T{AS) to B^B~ and B'^B^ are each 
50%. The branching fraction and /l are corrected for 
the slightly different reconstruction efficiencies in longi- 
tudinal and transversal polarizations. The statistical un- 
certainty on the signal yield is taken as the change in 
the central value when the quantity —2hiC increases by 
one unit from its minimum value. The significance is 
the square root of the difference between the value of 
— 21n£ (with systematic uncertainties included) for zero 
signal and the value at its minimum. In this calculation 
we have taken into account the fact that the floating /l 
parameter is not defined in the zero signal hypothesis. 

Figure [T] shows the projections of toes , , the ai 
invariant mass, J-' and Ti. for a subset of the data for which 
the ratio of the signal likelihood to the total likelihood 
(computed without using the variable plotted) exceeds a 
threshold that optimizes the sensitivity. 

A systematic uncertainty of 38 events on the signal 
yield due to the PDF parametrization is estimated by 
varying the signal PDF parameters within their uncer- 
tainties, obtained through comparison of MC and data 
in control samples. The uncertainty from the fit bias (7 
events) is taken as half the correction itself. Uncertainty 
from lack of knowledge of the ai meson parameters is 31 
events. We vary the SCF fractions by their uncertain- 
ties and estimate a systematic uncertainty of 12 events. 
A systematic uncertainty of 19 events from possible con- 
tamination by — > ai(1260)"'" a2(1320)~ background 
events is estimated with simulated MC experiments. The 
uncertainty due to cross feed between the signal and non- 
resonant backgrounds, evaluated with MC events, is 10 
events. Uncertainties of 1.4% and 3.6% are associated 
with the track efficiency and particle identification, re- 
spectively. Differences between data and simulation for 
the cos Ot variable lead to a systematic uncertainty of 




cos 9 

FIG. 1: Projections of (a) tties, (b) AE, (c) oi invariant 
mass (average of and m^- is shown), (d) T and (e) Ti. = 

I cos S I (average of jcosS^+l and |cosS^-| is shown). Points 

with error bars (statistical only) represent the data, the solid 
line the full fit function, and the dashed line the background 
component. These plots are made with a requirement on the 
signal likelihood that selects 25%-40% of the signal and 2%- 
5% of the background. 



2.5%. Assuming that 20% of oi decays proceed through 
the S-wave {tttt)^ channel [l8| . we estimate a systematic 
uncertainty of 6.8% from the difference in reconstruction 
efficiency between the P-wave (7r7r)p and S-wave com- 
ponents. The uncertainty in the total number of BB 
pairs in the data sample is 1.1%. The total systematic 
uncertainty, obtained by adding the individual terms in 
quadrature, is 12.9%. 

The main systematic uncertainties on /l arise from the 
fit bias (0.03), the variation of PDF parameters (0.08), 
the ai parametrization (0.04) and the non-resonant back- 
ground (0.02). 

In conclusion, we have measured the branch- 
ing fraction: S(B" ^ af ) x [B{at (37r)+)]2 = 
(11.8 ± 2.6 ± 1.6) X 10"^ and the fraction of longitudi- 
nal polarization = 0.31 ± 0.22 ± 0.10. Assuming 
that B{ai 7r^7r+7r+) is equal to B{af tt+tt^tt^), and 
that B{a+ (37r)+) is equal to 100% [S], we obtain 
B{B° ^ aj" ) = (47.3 ± 10.5 ± 6.3) x lO'^. The de- 
cay mode is seen with a significance of 5.0 a including 
systematic uncertainties. The measured branching frac- 



m ger 
in W, 



while they 



tion and longitudinal polarization are in general agree- 
ment with the theoretical expectations 
disfavor those in 2]. 
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